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Abstract 


A  computational  study  of  radiation  heat  transfer 
in  a  combustor  having  fuel  spray  combustion  Is  reported 
in  the  present  paper.  When  the  distribution  of 
radiatively  participating  species  Is  known  in  either  a 
normal  distribution  or  an  exponential  distribution  In 
any  direction  of  cylindrical  coordinates  with  respect 
to  a  referenced  location,  e.g.,  In  a  form.  F  «  f  exp  (- 
at-S;:-oz!),  where  f.  a,  b,  and  c  are  distribution 
constants  for  individual  species,  the  governing 
equation  or  radiation  heat  transfer  without  scattering 
in  such  a  system  can  be  solved  with  relative  ease  by 
using  the  herein  reported  method.  Several  results  from 
the  present  authors'  work  are  explained  that  are 
employed  for  Implementing  the  solution:  a  new 
coordinate  transformation  method  to  find  the  species 
distribution  along  individual  optical  paths  centered  at 
any  chosen  location  in  the  combustor;  a  new  formula  for 
the  adiabatic  flame  temperature  in  a  logarithmic 
function  expressed  In  terms  of  relevant  variables;  a 
new  integral  function  facilitating  solutions  for 
various  radiation  equations;  a  new  Inverse  error 
function  enabling  some  exact  solution  of  the  governing 
equation  of  radiation  heat  transfer,  etc. 

An  exact  solution  Is  found  for  both  optically  thin 
and  thick  media.  And  a  numerical  solution  Is  offered 
for  optically  medially  thin  medium.  The  present 
analytical  method  is  applied  to  radiation  heat  transfer 
In  a  direct  inject lon-type  diesel  engine  combustion 
chamber.  Some  representative  results  are  discussed 
that  have  been  obtained  from  a  parametric  analysis  of 
its  radiation  heat  transfer. 
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Plank  radiation  function 
Species  distribution  constant 
Plank's  constant 
Radiation  Intensity 
Boltzman  constant 
Radiation  heat  flux 
Optical  path 
Radius  of  cylinder 
Temperature 


•  ■.liability 


Components  in  cylindrical  coordinate 


Components  in  spherical  coordinate 

Volume  absorptance 

Wavelength 

Optical  depth 

Stefan-Boltzman  constant 


Subscripts 

a  Adiabatic 
b  Blackbody 
d  Detector 

1  1th  subrange  in  r-domaln 
J  Jth  subrange  In  A-domaln 
A  Spectral 

0  Foot  of  the  optical  path 


INTRODUCTION 


When  an  equation  of  radiative  heat  transfer  Is 
solved  for  predicting  the  heat  flux  Incident  upon 
various  locations  over  the  surface  of  a  combustor 
having  flame  plumes,  several  difficult  problems  are 
encountered.  In  order  to  briefly  review  the 
difficulties,  a  combustor  having  a  Jet  flame  Is 
considered  as  shown  in  Fig.  1.  The  heat  flux  through 
an  optical  path  (r,6,;)  centered  at  location  D  may  be 
calculated  by  solving  the  governing  equation  of 
radiation  neglecting  scattering  effect, 
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(1) 


F  -  f  exp(-ap-b4>*-czJ) 


(3) 


qe;  ’  ^  ^(0)  cose  dA 

0 

wnen  the  spectral  radiation  intensity,  I^(r)  Is 
represented  as 

rQ  <x(r’)  r' 

I ^ C r )  -  — - -  ebJ  exp(-  /  KJ(r")dr,')dr'  (2) 

r 

where  «r  is  the  spectral  volume  absorptanoe  at  (r,8,;) 
and  e  is  the  spectral  blackbody  radiation.  The  total 
radiation  flux  on  location  D,  then.  Is  calculated  by 
integrating  q  over  the  entire  hemispherical  volume 
faced  by  the  Ideation.  Among  the  main  difficulties  In 
implementing  the  above  computation  are  that  *.  can  be 
determined  only  when  the  distributions  of  radiatively 
participating  species  and  temperature  are  known  along 
the  individual  optical  paths  and  that  I.  Is  spectrum- 
dependent.  ' 


Nozzle,  t)(  p,$,z) 


Fig  I  A  Combustor  with  Jet  Flames. 


where  F  represents  either  fuel/air  ratio  or  soot 
concentration  or  CO,  or  H^O  concentration;  f,  a,  b  and 
c  are  constants  to  be  separately  determined  for  each 
medium  by  either  experimental  or  theoretical  means. 

The  expression  Is  considered  as  to  reasonably  well 
describe  the  species  distribution  of  non-axlsymmetrlc 
plumes  In  view  of  some  reported  results  from  both 
theoretical  and  experimental  studies  [6,7].  The  new 
expression  offers  various  advantages,  e.g.,  the 
expression  can  be  used  to  describe  the  spray  plume  In 
swirl  motion  by  including  Its  effect  In  the  in¬ 
direction  and  It  facilitates  the  sought  solution  of  Eq. 
(1)  [1].  When  the  species  distribution  with  respect  to 
a  referenced  position  In  the  combustor  Is  known,  the 
species  distribution  along  the  optical  path,  r,  of  any 

chosen  direction  (r,6,c)  and  location  D  (p.,*.,z.)  may 

d  dd 

be  found  as  [1], 

Fo  ’  fd  e*P(-C(r-r0)/ru]2  (*0 

where  f^,  r  and  r  are  functionally  related  to  6,  t, 

Pd*  *d‘  zd,0f'  a’  and  c'  Thls  dl3trlbutlon 
equation,  then,  may  be  used  for  determining  the  volume 

absorptance  according  to  the  reported  methods  C 3 , U , 5 ] 

as 

2 

-  <1  expt-c, (r-r1 )  )  (5) 

The  temperature  distribution,  T,  needed  In  the  above 
equation  and  for  Eq.  (1)  may  be  determined  by  using 
the  distribution  of  burned  fuel/alr  ratio,  w,  from  Eq. 
(3).  This  can  be  achieved  by  using  our  new  formula  for 
the  adiabatic  flame  temperature  [8], 

T  -  T  ♦  T.  lnu>  ♦  T,  Unu))2  (6) 

where  T  ,  T1  and  T2  are  functionally  related  to  the 
reaction  pressure,  the  Initial  mixture  temperature  and 
the  number  of  carbon  atoms  of  the  fuel  molecule.  Even 
with  the  above  rearrangement  and  simplification  of 
terms  in  Eq.  (1),  its  solution  was  not  readily  found 
until  some  new  mathematical  techniques  were  employed  as 
explained  next. 


In  order  to  obtain  a  better  solution  of  Eq.  (1), 
some  new  methods  were  proposed  for  overcoming  the 
abovement loned  difficulties  when  the  authors  computed 
radiation  heat  transfer  in  a  direct  Injection-type 
diesel  engine:  a  coordinate  transformation  technique  to 
find  the  species  distribution  along  Individual  optical 
paths  when  the  distribution  is  given  with  respect  to  a 
chosen  location,  e.g.,  the  Injection  nozzle  hole  [1];  a 
computation  method  of  finding  the  spectral  volume 
absorptances  of  soot  and  gas  for  their  mixtures  [1,2] 
by  using  reported  results  [ 3 ,  ** . 5 ]  and  new  analytical 
and  numerical  methods  of  solving  Eq.  (1),  which  are 
considered  in  the  present  paper.  This  paper, 
therefore,  may  partially  serve  as  a  summary  report  of 
the  authors'  modeling  study  of  radiation  heat  transfer 
In  a  combustor  having  flame  jets.  In  addition  to 
reporting  new  analytical  solutions  of  the  radiation 
neat  transfer  equation. 


ANALYTICAL  AND  NUMERICAL  SOLUTIONS 

The  solution  of  Eq.  (1)  may  be  considered  for  the 
following  three  cases:  (1)  the  optically  thin  medium; 
(2)  optically  thick  medium  and  (3)  optically  medially 
thick  medium.  The  solution  for  each  case  will  be 
explained  upon  further  simplifying  the  governing 
equation  as  follows.  Introducing  new  terms, 

kT  r 

K  -  <.*x  [3.9],  X  -  hc/kXT  and  a  -  f—  /  «  dr, 

o  a  nc  q  o 

where  h  is  Planck's  constant,  c  the  speed  of  light  and 
k  the  Boltzmann  constant,  Eq.  (1)  can  be  rewritten  for 
directional  radiation  Intensity  as, 

V  *  qec/c090' 


Elaborating  on  the  species  distribution  In  the 
combustion  chamber,  It  was  found  to  be  conveniently 
described  by  an  equation  in  cylindrical  coordinates 
with  respect  to  the  Injection  nozzle  [1]  when  the 
distribution  Is  known  In  either  a  norma  1  (or  a  skewed 
normal)  or  an  exponential  form  In  all  directions  of  the 
cylindrical  coordinate  system,  e.g.,  in  a  form, 
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During  the  process  of  simplifying  the  above  equation, 
the  authors  discovered  a  new  useful  integration 
function  [10], 


x.  -ax  m 
,  2  e  x 


/  -  dx  -  Fim.a.x^  -  F(m,a,x2), 

x  ex-1 


F(m,a,x)  • 


Z  !e-(n*a)x  I  [  L^ILL,  (n>a)'1-1xra-i]l. 
n-i  1*0 


where  r(m»1 )  -  m! 


Note  that  Eq.  (7)  can  also  be  used  for  determining 
a  solution  in  closed  form  of  the  black  body  radiation 


f  unction ,  F„A_  [  1 '  ] , 


IT  E  (T) 
bi 


*  -nx  ,  ,  2  , 

i  (x  *  Hr  ♦  H  *  • 

,  n  n  2  3 

n-1  n  n 


From  Eq.  (7),  since  one  can  find 


4  -ax 
x  e 


n-1  (n*a)' 


the  governing  equation  is  simplified  as 


i  360  ok  J_  ,'0  "o 

n5  0  (1.a/n)5 


With  the  definition  of  the  optical  depth,  t  -  /  dr. 


the  equation  may  be  rewritten  in  a  tidier  form, 


360  ok  ,  1  ,  o 

n-1  n  0 


Note  that  <  in  the  equation  is  assumed  constant  with 
respect  to  X,  which  appears  valid  when  soot  generated 
radiation  dominates  C'9].  The  solution  of  this 
equation  is  considered  for  the  previously  mentioned 
three  cases  as  follows. 


OPTICALLY  THIN  MEDIUM:  When  t  approaches  zero,  l.e., 
in  the  optically  thin  limit,  an  exact  solution  for  Eq. 
(61  was  found.  In  this  case,  Eq.  (8)  becomes 


.  360aw  ,  1  /  o  t5(1.  6kT, 

”  hc  n-1  n5  0  nhc 


For  finding  the  solution  of  the  above  equation,  since  T 
is  functionally  related  to  the  optical  depth,  t,  and 
the  burned  fuel/alr  ratio,  w,  and  since  u  may  be  found 


<3  3  'a)  *  u>,  e 


from  Eq.  (k),  Eq.  (6)  may  be  written  as 


T  •  T,  (in 
o  1 


r2)  ♦  T?  (In  u,  -  r2)2 


-  T3  -  T,  r  .  T2  r 


where,  T^  -  Tq  ♦  T,  In  w,  ♦  T?( In  u^)2  and 


T4  -  T,  ♦  2T2  tn(Ul). 


In  addition,  for  obtaining  the  solution  of  Eq.  (8),  it 
was  necessary  to  find  the  relationship  between  t  and  r. 
From  definition  and  Eq.  (5),  the  optical  depth,  t,  Is 
rewritten  as 


r  -c  (r-r.) 

t  -  <,  /  e  dr 

0 


“  2~  ^  c~  *erf  ♦  erf  (*'c“  r1 ) ) 


c 

erf  [kc7(r-r,)]  -  /  —  t  -  arf  (/c “  r,).  (11) 


In  view  of  Eq.  (10)  and  for  further  simplifying  Eq. 
Oi),  a  variable  transformation  was  employed  to  obtain 


“•  ^  (r-r, )  , 


i  — 7  /  —  [t  -  erf  (v£7  r  )] 

x,  X  2  °1  11 


d1  tT‘,lj  • 


Then  Eq.  (11)  may  be  written  as 


erf  (r)  -  t 


or,  r  -  erf  (t) 


where,  erf  represents  the  inverse  error  function, 
whose  exact  solution  is  found  as  follows  (the  details 
of  its  derivation  are  explained  in  the  Appendix): 


r  -  erf  1 (7) 


/»  x  —3  7**  —5  , 

r  (T  *  273T  T  *  ¥757  T  * 


Next,  the  above  equation  is  Introduced  into  Eq.  (10)  to 
find 


_  _  k  »  -3  7»*  -5  .2 

T  •  r3  -  IT  °  2717  T  *  ¥757  T  ••••) 


T2*a  _ w  -3  7»*  -5  . k 

2.3!  575 


*  —  2  M  w 2  _  .  — 

T3  ”  k  Tu  '  t_  kTTT  75  T2  1 


-  t3  -  T,  72  ♦  t6  7“ 


where,  T^  -  ^  Tu  and  Tfc  -  (-  ^  ^  T?) . 


From  Kg.  (lb),  one  may  also  find  the  following  In  oroer 
to  use  them  In  Eq.  (9): 


■  V%v 


W*1 1 


ii 


r  -  T,5  -  51,“  Tc  T2  *  (10T33  T52  *  5f3“  Tfi )  7*. 


‘3  3  5 


Tb  -  I3b  -  6T3  T5  x  ♦  (15T3“  T5"  *  6T3b  T6)  x\ 


Optically  thick  medium:  An  exact  solution  can  also  be 
found  In  an  optically  thick  medium:  assuming  the  volume 
absorptance  <  Is  0  fop  r  <  r.,  and  tor  rQ  <  r  < 
rQ*c,  one  can  rewrite  Eq.  ( 27  as 


Finally,  the  above  equations  and  Eq.  (12)  are 
Introduced  Into  Eq.  (9).  with  a  change  of  variables 

T  -  d^i-tj)  and  dT  -  1 /d1  d7,  to  obtain 
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‘e;  ‘  *  T  ' 

n-1  n  -d1i1 


,  5kx,T  6  5kTT  6  _ 

-  f(T  “  — * - — — )  —  - - —  .. -  f 

d1  3  nhc  nhcd1 


r  »e  r’ 
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This  equation  Is  further  solved  to  find  the  result, 
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Optically  medially  thin  medium:  When  the  system  holds 
an  optically  medially  thin  (or  thick)  medium,  Eq.  (2) 
cannot  analytically  be  solved  as  the  above  discussed 
cases;  rather,  a  numerical  method  was  employed  for  It. 
Without  going  Into  the  details  of  Its  derivation,  but 
by  knowing  that  the  volume  absorptance  <(r)  can  be 
found  from  Eq.  (5),  the  following  result  was  used  for 
the  solution. 


-^<’5T3T52*6T35  V7^7  (16> 

Consequently,  Eq.  (16)  can  be  used  for  finding  an  exact 
solution  of  the  equation  of  radiation  for  an  optic  ’ ly 
thin  medium  In  a  combustor  with  flame  plumes  as. 
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*here»  *c  is  the  volume  absorptance  defined  In  the 
1th  subrange  of  r-domaln  and  the  Jth  subrange  of  X- 
domain.  Then, 
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where,  dx  -  v  dr.  The  numerical  results  of  Eq. 

(  '  9 )  were'obtalfteo  using  Causslan  Integration  method. 


RESULTS  AND  CONCLUSION 

The  present  computational  method  enables  an 
extensive  parametric  study  of  radiation  heat  transfer 
in  a  combustor  with  flame  Jets  [2J.  Since  this  paper 
is  to  primarily  report  the  analytical  methods  of 
computation  by  using  the  authors’  several  new  findings 
reported  elsewhere  [1,  2,  5,  '0,  1 1  j .  only  some 
represent  at  1 ve  results  obtained  for  a  direct  Injection- 
type  liesel  engine  are  presented  here. 


Note  that  since  our  computer  program  to 
numerically  Implement  Eq.  (1)  was  developed  prior  to 
finding  the  exact  solutions  that  are  reported  here  and 
since  some  comparison  between  the  results  obtained 
using  the  program  and  those  using  the  exact  solutions 
was  satisfactory,  the  results  discussed  next  were 
mainly  obtained  using  our  available  computer  program. 

In  the  computation,  in  order  to  be  consistent  and  to 
have  better  mutual  comparisons,  the  same  computational 
conditions  employed  for  the  previous  results  are 
considered  here.  The  numerical  values  of  f  in  Eq..(3) 
were  for  soot,  HO,  CO,  and  fuel/alr  ratio  8  x  10  , 

0.01,  O.Oi  and  170,  respectively:  the  distribution 
constants,  a,  b,  and  c  were,  for  these  species,  0.6, 
2.36  and  2.0,  accordingly.  Note  again  that  these 
constants  and  the  distribution  made  in  each  coordinate 
of  Eq.  (3)  can  be  determined  by  computational  or 
experimental  methods,  or  arbitrarily  chosen  for 
parametric  studies.  The  engine  details  for  the 
comgutatlon  were:  the  chamber  surface  temperature, 

500  K;  the  number  of  spray  plumes,  H ;  the  piston 
radius,  R  -  U.92  cm;  the  combustion  chamber  bowl 
radius,  3. MM  cm;  the  compression  ratio,  23:  the  engine 
speed,  1000  rpm;  the  surface  spectral  emisslvlty,  0.95: 
the  fuel  composition,  the  time  of  computation 

in  a  cycle,  M0  degrees  arter  top  dead  center.  The 
results  are  shown  in  normalized  forms:  normalized 
radiation  heat  transfer,  Q/oT  ,  where  T  is  2M00°K  and  o 
is  the  Stefan-Boltzmann  constant:  normalized  direc¬ 
tional  intensity,  I  /xaT  :  normalized  spectral 
radiation  heat  transfer,  Q^/aT  . 


I  is  7  Directional  Radiation  Heat  Flux  Incident  on  Cylinder 
Head  Surface  in  Radial  Direction 


The  radiation  heat  transfer  through  a  solid  angle 
incident  on  tne  cylinder  head,  are  shown  in  Fig.  2. 

•hen  tnere  are  four  flame  Jets  in  the  combustion 
chamber,  radiation  heat  transfer  flux  incident  on  the 
locations  along  the  radial  direction  varies  depending 
jpon  the  direction  in  the  hemispherical  volume  of 
integration.  The  schematic  drawing  in  the  figure  shows 
t.ne  spherical  coordinates  identifying  the  glrectlon  gf 
solid  angles.  Note  that,  for  cases  8-60  and  76.7  , 
neat  transfer  on  r/R  <0.1  is  largely  from  plume  #1  and 
that  beyond  r/R  -  0.1  from  both  plumes  #1  and  #2.  The 
directional  radiation  intensity,  I  was  combined  with 
ff  spectral  emisslvlty  of  the  chamber  wall,  Cj  -  0.95, 


Fig.  3  Radiation  Heat  Transfer  through  Cylindrical  Head. 


to  obtain  radiation  heat  transfer  through  the  cylinder 
head  by  integrating  as,  Q  -  /  I  cose  sin  8  d6  dc 
(Fig.  3).  The  net  heat  transfer^computed  with 
inclusion  of  wall  emission  at  500K,  reflection  from  the 
opposite-side  walls  are  represented  by  "Net";  those 
obtelned  excluding  the  reflection  are  denoted  by  "No 
reflect  ion" :  and  the  heat  flux  incident  on  the  surface 
without  including  the  wall  effect  ere  shown  by  "Spray 
plune".  Explaining  the  results,  the  rapid  reduction  In 
heat  transfer  around  r/R  -  0.7  Is  caused  by  the 
presence  of  the  piston  bowl.  The  double  hump  In  the 
plot  may  be  surprising.  In  order  to  clarify  this,  an 
additional  computation  made  for  solid  angle  of  zenith 
angle,  8  -  39.7°,  at  various  azimuthal  angles  and 
locations  on  the  cylinder  head,  is  shown  In  Fig.  M. 

From  the  results,  one  finds  that  while  the  radiation 
heat  flux  Incident  on  locations  around  the  nozzle  (r/R 
-  0)  Is  from  all  of  the  four  plumes,  the  flux  on  those 
near  r/R  -  0,1  is  from  one  plume,  l.e.,  the  nearest 


Fig  4  Directional  Intensity  of  Radiation  Incident  on  Cylinder  Head 
through  Various  Directions  of  Solid  Angle. 
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plume  to  the  location.  If  the  computation  is  made  for 
a  single  plume,  the  heat  transfer  should  have 
asymptotically  been  approaching  zero  at  r/R  -  0.  This 
trend  is  qualitatively  compared  with  the  experimental 
results  obtained  by  using  a  single  plume  [12]. 

Since  the  distribution  constants  being  used  for 
Eq.  (3)  in  the  present  analysis  were  chosen  based  on  an 
in-cylinder  soot  measurement  [13].  it  seems  of  interest 
to  compute  radiative  heat  transfer  contritJted  by  in¬ 
cylinder  soot  and  gaseous  species,  l.e.,  mainly  CO^  and 
HjO.  Fig.  5  shows  their  spectrum-resolved  radiation 


Normalltrd  Number. 


fig  5  Radiation  Meat  Transfer  Contributed  by  Soot  and  Gaseous 
Species  and  their  Combination. 

heat  transfer  smoothed^over  individual  wave  bands  In 
normalized  form,  Q^/oT  .  with  respect  to  normalized 
wave  number,  6  um/5.  It  appears  clear  that  thermal 
radiation  transmission  caused  by  the  presence  of  soot 
almost  dominates  the  entire  radiation  process  and  the 
radiation  by  gaseous  species  is  very  small.  The  great 
portion  of  radiation  by  gas  found  around  6um/X  •  1.8  is 
explained  by  the  strong  emission  bands  of  H?0  at  A  * 
2.7um  and  CO,  at  1  -  2.7  and  9.3  urn.  Since  the 
spectral  ab30rptance  of  a  surface  can  be  Incorporated 
with  results  like  those  in  Fig.  5,  the  analysis  may  be 
made  for  combustors  having  various  surface  coatings 
with  Known  spectral  eml33lvity. 

The  main  issue  associated  with  the  use  of  the 
present  analytical  method  of  radiation  heat  transfer 
may  be  Justification  of  the  species  distribution 
described  by  Eq.  (3)-  Since  the  distribution  widely 
varies  depending  upon  the  fuel  injector  and  combustor 
condition,  It  Is  difficult  to  exactly  describe  the 
variations  by  a  single  equation.  Reported  results 
suggest,  however,  the  individual  species  distribution 
may  be  simplified  by  using  an  exponential  or  normal  (or 
skewed  normal)  distribution  in  each  coordinate 
component  of  a  cylindrical  system,  e.g.,  Eq.  (3),  It 
Is  found  In  the  present  study  that  the  radiation  heat 
analysis  is  greatly  simplified  when  the  distribution  is 
made  In  such  a  form  without  resorting  to  the  rather 
complex  zonal  method  or  the  use  of  the  geometric 
factors.  Further,  the  use  of  the  present  method 
readily  enables  a  parametric  analysis  of  radiation  heat 
transfer,  e.g.,  the  selection  of  an  Injector  suitable 
for  a  chosen  combustor  system. 


The  present  study  has  been  sponsored  by  the  U.S. 
Army  Research  Office,  Contract  No.  DAAG  29  -83-K-0093 
(Scientific  Program  Officer,  Dr.  David  M.  Mann). 

REFERENCES 

1.  Chang,  S.  L.  and  Rhee,  K.  T. ,  "Coordinate 
Transformation  Method  for  Radiation  Heat  Transfer 
Prediction  in  Soot  Laden  Combustion  Products," 
being  reviewed  for  publication. 

2.  Chang,  S.  L. ,  Yang,  X.  L.  and  Rhee,  K.  T.,  "A 
Parametric  Analysis  of  Radiation  Heat  Transfer  in 
Direct  Injection  Diesel  Combustion,"  International 
Symposium  on  Diagnostics  and  Modeling  in  Reciprocating 
Engines,  September  U— 6,  1985,  Tokyo,  Japan. 

3.  Dalzell,  W.  H.  and  Sarofim,  A.  F. ,  "Optical 
Constants  of  Soot  and  Their  Application  to  Heat 
Flux  Calculation,”  J.  of  Heat  Transfer,  Trans,  of 
ASME,  Vol .  91.  pp.  100-109,  1969. 

9 .  Lee,  S.  C.  and  Tien,  C.  L. ,  "Optical  Constants  of 
Soot  in  Hydrocarbon  Flame,"  18th  Symposium  (Int'l)  on 
Combustion,  The  Combustion  Institute,  pp.  1159-1165, 

'981  . 

5.  Edwara,  D.  K.  and  Balakrlshnan,  A.,  "Thermal 
Radiation  by  Combustion  Gases,"  Int’l  J.  of  Heat 
Mass  Transfer,  Vol.  16,  p.  25-90,  1973. 

6.  Takeuchi,  K.,  Senda,  J.  and  Shlkuja,  M.  ,  "Transient 
characteristics  of  Fuel  Atomization  and  Droplet 

Size  Distribution  in  Diesel  Fuel  Spray,"  SAE  Paper 
830999,  1983. 

7.  Kuo,  T.  W.,  Yu,  R.  C.  and  Shahed,  S.  M.,  "A 
Numerical  Study  of  the  Transient  Evaporating  Spray 
Mixing  Process  in  the  Diesel  Environment,"  SAE 
Paper-831735,  1983. 

8.  Rhee,  K.  T.  and  Chang,  S.  L. ,  "Empirical  Equations 
for  Adiabatic  Flame  Temperatures  for  Some  Fuel-Air 
Combustion  Systems,"  Combustion  Science  and 
Technology,  Vol.  99,  pp.  75-88,  1985. 

9.  Grosshandler,  W.  L.  and  Modak,  A.,  "Radiation  from 
Nonhomogeneous  Combustion  Products,"  18th  Symposium 
(Int’l)  on  Combustion,  The  Combustion  Institute, 

pp.  601-609,  1981. 

10.  Chang,  S.  L.  and  Rhee,  K.  T, ,  "A  Useful  Integral 
Function  and  Its  Application  in  Thermal  Radiation 
Calculations,"  International  Communication  of  Heat 
and  Mass  Transfer,  Vol.  10,  No.  9,  pp.  329-333.  1983. 

11.  Chang,  S.  L.  and  Rhee,  K.  T. ,  "Blackbody  Radiation 
Functions,"  International  Communication  of  Heat  and 
Mass  Transfer,  Vol.  11,  pp.  951-955,  1989. 

12.  Tsubol,  T.,  Fukushlma,  M.,  Sato,  T.  and  Ogurl,  T. , 
"Thermal  Radiation  During  Spray  Combustion  Behind 
Reflected  Shock  Waves,"  International  Symposium  on 
Diagnostics  and  Modeling  In  Reciprocating  Engines, 
September  9-6,  1985,  Tokyo,  Japan. 

13.  Klttleson,  D.  B. ,  Du,  C.  J.  and  Pipho,  M.,  "Total 
Cylinder  Sampling  from  a  Diesel  Engine:  Part  II  - 
Particulate  Measurement",  SAE  Paper  830293,  1987. 

'9.  Papoulls,  A.,  Probability,  Random  Variables,  and 


6 


Stochastic  Processes,  McCraw  Hill  Book  Co.,  New  x"(v)  -  —  fx'(y)]  -  —  [x' (y)  X^-]”' 

York,  1 965 .  dy  dx  dx 

15.  Schneider,  P.  J.,  Conduction  Heat  Transfer, 

Addison-Wesley  Publishing  Co.  Inc.,  Reading,  Mass.,  x  2x2 

1 95c  _  *  j  x  e  , 


16.  Abramouitz,  M.  and  Stegun,  I.  A.,  Handbook  of 
Mathematical  Functions,  National  Bureau  of 
Standards,  Applied  Mathematic  Series-55,  U.  S. 
Government  Printing  Office,  Washington,  D.C.,  196k. 


Appendix : 

AN  INVERSE  ERROR  FUNCT I ON 


The  error  function,  erf(x),  although  defined  in 
several  ways,  is  essentially  equivalent  to  the 
following  expression: 


erf ( x ) 


/x  0 


_2 

/x 
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n-0 


(-1)" 


2n*1 


_x _ 

n !  ( 2n ♦  1 ) 


(A-1) 


The  equation  is  the  integral  of  the  so-called  Gaussian 
or  normal  function  and  occurs  frequently  in  the  study 
of  the  general  theory  of  probability  [Ik],  the  analysis 
of  transient  heat  flow  in  a  seml-lnflnlte  solid  [15], 
and  the  computation  of  radiation  heat  transfer  with 
radiatively  participating  media  [10],  etc.  For  meeting 
the  needs  of  such  studies,  the  error  function,  its 
derivatives  and  Integrals  have  been  tabulated  [16].  As 
found  in  the  text,  it  was  needed  to  obtain  the  Inverse 
form  of  tne  error  function,  i.e.,  in  Eq.  (A-1),  the 
value  of  x  for  a  corresponding  erf(x).  The  following 
derivation  shows  a  new  inverse  error  function  in  closed 
form: 
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Consequently,  one  can  rewrite  Eq.  (k)  as 


When  the  error  function,  erf(x),  is  3et  equal  to 
y,  i.e.,  y  -  erf(x),  its  inverse  form  may  be  expressed, 
for  the  convenience  of  discussion, 

x  -  erf-1 (y )  (A-2) 


x  -  fre(y)  -  ~  (y  ♦  — j-  y3 
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or , 


x  -  fre(y).  ( A-3> 

By  Taylor's  expansion  method,  the  Inverse  error 
function  can  be  expressed  as 
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where,  x(n)  (0)  is  the  nth  derivative  of  x 
to  y  at  y  -  0.  In  order  to  find  xCn'  (0), 
following  steps  are  taken. 
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A  close  Investigation  of  the  above  results  Indicates 

that  a  term  ‘S’  is  repeated  in  each  derivative  and 

that  the  even  derivatives  become  zero  at  y  -  0.  This 
leads  to  writing  Eq.  (5),  the  Inverse  error  function, 
in  a  tidier  form. 


fre(y)  -  E  f  (0) y)2"*1  /(2n+1 )  I  |y|  <  1  (A-6) 

n-0  n  2 

where,  fn(0)  is  obtained  from  the  reciprocally  defined 
function,  fn(t),  expressed  as 

f  ( t)  -  1 ,  and 
0 

f  . . ( t )  -  f ' ' ( t )  ♦  2(kn  ♦  3)  X  f'(t) 
n*  i  n  n 

♦  2[2n*1  ♦  k( 2n*1 )  (n*l )  x2]  f  (t). 

n 

(n  -  0,1,2,...) 

Among  the  analytical  properties  of  the  new  Inverse 
error  function  are 
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